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ABSTRACT 

We present an analytic model for the sizes of voids in the galaxy distribution. Peebles and 
others have recently emphasized the possibility that the observed characteristics of voids may 
point to a problem in galaxy formation models, but testing these claims has been difficult with- 
out any clear predictions for their properties. In order to address such questions, we build a 
model to describe the distribution of galaxy underdensities. Our model is based on the "excur- 
sion set formalism," the same technique used to predict the dark matter halo mass function. We 
find that, because of bias, galaxy voids are typically significantly larger than dark matter voids 
and should fill most of the universe. We show that voids selected from catalogs of luminous 
galaxies should be larger than those selected from faint galaxies: the characteristic radii range 
from ~ 5-10 hr 1 Mpc for galaxies with absolute 7--band magnitudes M r — 5 log h < — 16 to 
—20. These are reasonably close to, though somewhat smaller than, the observed sizes. The 
discrepancy may result from the void selection algorithm or from their internal structure. We 
also compute the halo populations inside voids. We expect small haloes (M < 10 11 M ) to be 
up to a factor of two less underdense than the haloes of normal galaxies. Within large voids, 
the mass function is nearly independent of the size of the underdensity, but finite-size effects 
play a significant role in small voids ( < 7 h~ x Mpc). 



Key words: cosmology: theory - 
functions 



large-scale structure of the universe - galaxies: luminosity 



1 INTRODUCTION 

One of the key predictions of any cosmological structure formation 
model is the distribution of matter on large scales. We now know 
that the cold dark matter (CDM) paradigm can account for many 
of the observed characteristics of the galaxy distribution. In this 
picture, the universe contained tiny (gaussian) density fluctuations 
at the time the cosmic microwave background was last scattered. 
Bound structures assembled themselves through gravitational in- 
stability around these perturbations in the relatively recent past. In 
the past two decades, with the advent of high-resolution numerical 
simulations, it has become possible to follow this picture through 
the formation of massive galaxies and clusters. The CDM model 
accurately describes the abundance and clustering of collapsed ob- 
jects from dwarf galaxies to rich galaxy clusters over a wide range 
of redshifts (although the details of galaxy formation itself remain 
somewhat mysterious) as well as the distribution of neutral gas in 
the intergalactic medium. 

These systems all correspond to density peaks in the mat- 
ter distribution (with the exception of the lowest column den- 
sity Lya forest absorbers). For a variety of reasons, the other 
end of the density distribution - underdense voids - has received 
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considerabl y less attention, despi t e their long o bservational his- 
tory iGregory & Thompson 1978; Kirshner et alJll98 If) and their 
place as the most visually striking features of the galaxy dis- 
tribution. This is largely because voids subtend enormous vol- 
umes and so require large surveys to garner representative sam- 



pies. Although voids have been found 


in every redshift 


survey 


Ide LapDarent et al. 1986; Vogelev et al. 


11994 lEl- Ad et alJll997l 


El- Ad & Piran 2000; Miiller et al. 2000; Hovle & Vogelev 2002) 



the first statistically significant sample came only with the 2dF red- 
shift survey <Hovle & VogelevlEool) . Now, with the DEEP2 red- 
shift survey and the Sloan Digital Sky Survey, it is even possible to 
constrain the evolution of voids over the redshift interval z ~ 1-0 

IConrov et all2005h. 

iHovle & Vogelevl |2004) presented the most complete search 
for voids to date. They found that voids with characteristic radii 
R « 15 h^ 1 Mpc fill ~ 35% of the universe. But their search il- 
lustrates a second difficulty in studying vo ids: how to define and 
identify them precisely (and meaning fully). lEl- Ad &Pirad 1 19971) 
proposed the Voidfinder algorithm based on separating the ob- 
served galaxies into "void" and "wall" populations and building 
voids around gaps in the wall population (see also Hovle & Vogelev 
l2002t) . While clearly defined for any given observational sam- 
ple, the results can nevertheless be difficult to interpret in rela- 
tion to the underlying physical quantities of interest. For exam- 
ple, the distribution of observed sizes depends on the galaxy sam- 
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pie (intrinsically brighter galax ies yield large r voids) as well as 
the search algorithm I Hovle & Vogelev 2004 restrict their search 
to radii greater than 10 h^ 1 Mpc, for example). Smaller voids are 
difficult to pinpoint because of confusion with random fluctua- 
tions in the galaxy distribution. We will nevertheless follow this 
approach and define a void to be any coherent region where the 
galaxy density falls below some threshold. Note that this differs 
from many other studies (e.g..lEinasto et alJll989tlGottl6ber et alj 
l2003t ISheth& van de Wevgaert! 120041) that require a void to be 
completely empty of galaxies. These voids, which are much smaller 
than the voids we consider here, do not correspond to the voids de- 
tected by eye in the galaxy distribution. 

The void phenomenon is also relatively difficult to study the- 
oretically. The simplest model is a spherical tophat underdensity. 
The early evolution of such a system is well-described by spher- 
ical expan sion (the a nalog of the well-known spherical collapse 
model: |Peebiesl ll980). Underdensities expand in comoving units, 
gradually deepening, until they reach "shell-crossing," when the 
center is evacuated and individual mass shells cross paths. At this 
stage, the initial spherical expansion model breaks down. Later, the 
voids continue to e xpand rel atively slowly in a sel f-similar fashion 
iFillmore & Goldreichll98llSuto et alll984lBertschingeJl985h . 
Unfortunately, these dark matter models are idealizations in that 
real voids are observed only through the galaxy distribution. Be- 
cause galaxies are biased relative to the dark matter, we must ask 
what kind of physical systems observed voids actually represent. 
They are in fact nearly empty of galaxies - but does that require 
shell-crossing, or can they be at an earlier evolutionary stage? 

Another difficulty is that the large size of voids restricts 
the usefulness of numerical simulations. Early efforts focused 
on understanding the dynamics of indivi dual voids, which did 
not require particularly high resolu tion iDub inskiet"aTlll993t 
van de Wevgaert & van KampenlfT99l . But placing voids in their 
proper cosmological context demands both large volumes and high 
mass resolution - the latter because we mus t resolve t he galaxies 
from whose absence we identify voids (e.g.. lGoldberg & Vogelev! 
2004). Only recently have iV-body simulations of the required 
dynamic range become practical. This has allow ed the first sys- 
temat ic studies of the structure of large voids l Gottlo ber et alJ 
120031) as well as of simulate d voids in the dark matter dis- 
tribution IColberg et alj 120051) and in t he galaxy distribution 
(Mathis & White 2002;: lBenson et all2o6H) . Nevertheless, there are 
still no full hydrodynamic simulations of the void phenomenon: 



galaxy properties are currently 


' determined througl 


models (Mathis & Whitel2002t 


iBenson et alJ2003f) 



There are a number of reasons to study the void phenomenon. 
Early attention focused on usin g the observed vo ids t o con strain 
cosmo logical parameters. Blum enthal et alJll992t) and lPiran et alJ 
( 1993) argued that the scale of typical voids would depend on the 
matter power spectrum, just as the abundance of clusters does. It 
was quickly realized that the void distribution seemed to have much 
more large-scale power than collapsed objects - most obviously, 
although voids with sizes > 15/t -1 Mpc are not uncommon, col- 
lapsed objects do not reach the same mass scales. This is puzzling 
because an underdense region reaches shell-cros sing only a fter the 
equivalent overdensity would virialize. lFriedmann & Pirarl 1200 1!) 
pointed out that a solution might lie in a proper treatment of the 
galaxies used to define voids: galaxy bias could lower the required 
dark matter underdensity, allowing voids to be larger for a given 
amount of CDM power. They argued that the observed population 
of large voids preferred a ACDM model. 

Recently, interest in voids has focused on their role in galaxy 



formation. IPeebleJ 1200 lh argued that (in the CDM model) voids 
should be populated by small dark matter haloes but that the ob- 
served voids appear to lack faint galaxies as well as bright ones. Is 
this a fundamental problem for the CDM paradigm, or does it indi- 
cate that galaxy formation proceeds differently in voids? One pop- 
ular explanation is that photoheating during reionization may have 
suppressed the formation of dwarf galaxies in low-density environ- 
ments iTullv et alJl2002UBarkana & Loebll2004ICerlEoo3) . help- 
ing to clear faint galaxies from the voids. To study these problems, 
recent attention has focused on the variation of the galaxy luminos- 
ity function with the la rge-scale environment. ICroton et al J 1 2005 ) 
and lHovle et alj 120051) (see also lGoldberg et alj|2005h found that 
the characteristic galaxy luminosity and the galaxy density both de- 
crease significantly in low-density environments but that the faint- 
end slope remains nearly constant. The latter may be surprising 
given the expected variation in the halo mass function with envi- 
ronment. 

Before we can answer any of these questions, however, 
we require a proper u nder standing of voids i n CDM models. 
iMathis & White! 120021) and lBenson et alJ 120031) made important 
strides forward with their studies of voids in A-body simulations. 
They claimed that current semi-analytic galaxy formation models 
predict void properties si milar to those ob served and argued that 
the conflict pointed out by Peebles 1 200 1 ) was illusory. But such 
studies are still limited by their dynamic range and by the (nec- 
essarily complex) galaxy formation prescriptions imposed on the 
dark matter haloes. 

The goal of this paper is to produce a straightforward an- 
alytic model of the void distribution and of galaxy populations 
within voids. Such a model will sharpen our understanding of 
voids in a CDM model and generate a baseline prediction with 
which we can contrast their obse rved properties. We will build on 
ISheth & van de Wevgaert! J2004), the most compelling theoretical 
model of voids to date. They used the excursion set formalism, 
which reproduces the abundance of collapsed haloes extremely 
well, to predict the distribution of void sizes. However, they re- 
quired voids to reach shell-crossing and defined their properties in 
terms of the dark-matter underdensity. As a result, they predicted 
characteristic void radii much smaller than those observed. Our 
main goal is to modify their approach so as to describe voids in the 
galaxy distribution. Along the way, we will also be able to predict 
the halo populations within voids and quantify the claimed discrep- 
ancy with observational results. 

The remainder of this paper is organized as follows. In [J2] we 
examine the nonlinear evolution of the void density. Then, in f^l we 
show how to compute the linearized underdensity of voids with a 
specified galaxy underdensity, and we briefly discuss the expected 
galaxy populations inside voids. In Sj4] we show how to compute 
the cosmological abundance of large galaxy voids and compare our 
results to observations. Finally, we conclude in ^5| 

In our calculations, we assume a cosmology with Q m =0.3, 
n A = 0.7, Q b = 0.046, H = 100/1 km s" 1 Mpc" 1 (with 
h — 0.7), n = 1, and ag = 0.9, consistent with the most recent 
measurements (Sper gel et all2003h . 



2 THE GRAVITATIONAL EXPANSION OF VOIDS 

Because the principal quantity of interest is the physical volume 
of voids, we first describe how they expand beyond their initial 
comov ing size. We will follow the method of Friedm ann & Pirarl 
(2001). We consider the expansion of a tophat density perturba- 
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tion that begins with a (negative) density perturbation Si inside a 
physical radius Ri at an initial time U (corresponding to zi). If we 
assume zero peculiar velocity perturbation and consider only the 
epoch before shell-crossing (so that the mass inside of each shell 
is conserved), conservation of energy yields the equation of motion 
JPeebleJl980l) 



R 2 (t) = Hl[—Q.o{l + R~5i 

+fio(l + Zif{l + 8i)R 3 /R(t) + n A R{t) 2 



(1) 



where the first term on the right hand side is (twice) the initial to- 
tal energy of the shell, the second term is (twice) the gravitational 
potential energy at time t, and the last term is (twice) the excess 
energy from the cosmological constant. We solve this equation for 
the physical radius of the shell as a function of time or redshift. 
The solution is such that the void expands in comoving units once 
nonlinear effects set in. We define 

R(l + z) 



'I 



(2) 



Ri(l+Zi) 

to be the ratio of the comoving size at redshift z to its initial co- 
moving size. The real fractional underdensity due to gravitational 
expansion is thus 1 + 8 — 77 ~ 3 . Note that this is independent of Ri. 

Equation Q includes nonlinear expansion, but most of the 
following will be phrased in terms of the equivalent linear den- 
sity 8 L extrapolated to the present day. (For clarity, we will denote 
all linear-extrapolated densities with a superscript "L" in the fol- 
lowing.) In order to transform to physical densities and scales, we 
require ri(S L ) or equivalently 8(8 L ). We therefore also compute 
the density as if linear theory were always accurate. Any pertur- 
bation can be divided into growing and decaying modes; in linear 
theory the growing mode (which is of course th e only component 
surviving to the present day) obeys iHeathl 1977b 

r _ 3SiH(z 



du 



u + 1 
H 3 (l/u- 



(3) 



where H 2 {z) — S7o(l + z) 3 + Qa in a flat universe. The con- 
stant 3/5 comes from our choice of initial conditions: a perturba- 
tion with zero peculiar velocity has 60% of its amplitude in the 
growi ng mode and the remainder in the decaying mode IPeeblea 

We show the relation between the linearized and physical un- 
derdensities with the solid line in Figure Q At first, while linear 
evolution is accurate, 8 « 8 L , but then the physical underdensity 
flattens out and even a small deepening of the void corresponds 
to a large increase in \5 L \. The dashed line in Figure \\\ shows a 
fitting function for S(S L ) from eq. (18) of Mo & White 1 19961) . 
Although those authors were interested in 8 > 0, the formula 
is accurate to a few per cent in the underdense regime as well. 
For reference, shell-crossing occurs at 8 L = —2.8, correspond- 
ing to a physical underdensity 8 = —0.81. Beyond this point, our 
model for r\ breaks down because the enclosed mass no longer re- 
mains constant. The expansion then proceeds self-similarly, with 
the comoving radius R oc a 1 ^ 3 in an Einstein-de Sitter univers e 
iFillmore & Goldreicrll98llSuto et all 1 984 iBertschingeil 19851) . 
significantly slower than the expansion before shell-crossing. As a 
result, the void volume fraction will be dominated by regions that 
have not yet (or have just) reached shell-crossing 1 Blumenthal et al. 
Il992l) : deeper perturbations are rarer in the CDM model and can be 
neglected because they expand so slowly. We therefore do not con- 
sider the later evolution in any more detail. 

For completeness, we also note that Figure Q is nearly inde- 
pendent of the input cosmological parameters; differences in the 



-0.4 




<5 L 



Figure 1. Physical dark matter underdensity <5 as a function of the linearized 
underdensity 8 (solid line). Note that our model actually breaks down for 
S < —0.81, when shell-crossing occurs (marked by the horizontal dotted 
line). The dashed line shows the fitting function o flMo& White! 119961) . 



cosmology manifest t hemselves in the growth rate of 8 L rather than 
in the function 8(S L ) iFriedmann & Piran 2001). 



3 DEFINING VOIDS 

It may now seem straightforward to compute the void distribution: 
one simply chooses a physical density threshold 8, transforms to 
8 L , and uses the excursion set formalism to determine the mass 
distribution of objects with the specified linearized underdensity 
ISheth & van de Wevgaertll2004l) . While such a procedure is per- 
fectly well-defined, it has one crucial shortcoming: in observed 
samples, we do not define voids in terms of their dark matter den- 
sity but in terms of their galaxy density. Thus we actually want to 
relate the physical galaxy underdensity <5 ga i within some region to 
8 L . That is the task of this section. 



3.1 The galaxy mass function 

We will use the exc ursion set formal ism to derive the dark matter 
halo mass function I B ond et alJ 119911) . In the simplest such model, 
a halo of mass forms whenever the smoothed density field first 
exceeds a linear-extrapolated density 8^ that is fixed by the physics 
of spherical halo collapse (8^ — 1.69 for objects collapsing at the 
present day; Peebles 1980). To solve such a problem, we consider 
diffusion in the (8 L , a 2 ) plane, where a 2 is the variance of the 
density field smoothed on scale m^. Trajectories begin with S L — 
at a 2 =0 (i.e., they must have the mean density on infinitely large 
scales). They then diffuse away from the origin as density modes 
are added on smaller scales. The problem is simply to compute the 
distribution of a 2 (or equivalently mh) at which these trajectories 
cross the absorbing barrier 8^ for the first time. From this first- 
crossing distribution, the mea n comovin g number density of haloes 
with m asses mj, ± dm^ /2 is IPress & SchechteJl974tlBond et all 
Il99ll) 
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Table 1. Mass Thresholds and Survey Depths 



M r - 5 log h 


M b j - 5 log h 


n(10- 3 h 3 Mpc~ 3 ) 


m min (PS; 10 11 Mq) 


m min (ST; 10 11 M ) 


-16 


-16.3 


58.6 


1.0 


0.73 


-18 


-17.8 


27.9 


2.5 


1.7 


-19 


-18.7 


15.2 


5.2 


3.4 


-20 


-19.5 


5.89 


15 


10 


-21 


-20.4 


1.11 


90 


60 



Number densities are computed from the Blanton et al. i 2003) r-band and lCroton et all 120051) fej-b and luminosity 
functions. Columns 4 and 5 give the minimum halo mass required to match these densities, with the Kravtsov et al. 
|2004|) {N(m h )), for the Press-Schechter and Sheth-Tormen mass functions, respectively. 



n h (m h ) = \ 



2 p 



d In a 



dlnm^ 



si 

— exp 
a 



(Si 



2<r 2 



(4) 



where p is the mean matter density. In a region with linearized un- 
derdensity 8^ and mass M v [with M v = (4tt/3 )p(R v /v) 3 for a 
void], the comoving number density of haloes is (Bond et all 199lt 
lLacev & Colell993}) 



n h (m h \8 v , M v ) 



d In a 



x exp 



d In 

~2(a 2 -a 2 ) 



— 8y) 

(a 2 - a 2 ) 3 / 2 



(5) 



where a v = cr(M v ) and of necessity a > <j v . This conditional 
mass function follows from an identical diffusion problem, but in 
this case the trajectories begin at (8%, a 2 ,) rather than the origin. 

Although equation J4} provides a reasonable match to halo 
abundances in cosmological simulations, it is by no means perfect. 
Shet h & TormerJ 1 19991) and ljenkins eTai] 1200 ll) provide more ac- 
curate fits to the simulation results. We will nevertheless use the 
standard Press-Schechter abundances in the following. The main 
reason is that, although the Sheth-Tormen mass function can also be 
motivated by a diffusion problem with the absorbing barrier fixed 
by ellipsoidal collapse, there is no analytic form for the conditional 
mass function because the barrier ch anges i ts effective shape with a 
shift of the origin. Isheth & Tormerl j2002) showed that the condi- 
tional mass function could be approximated through a Taylor series, 
but in that case the corresponding unconditional mass function dif- 
fers slightly from the usual Sheth-Tormen fo rm (which pro vides the 
best match to numerical simulations). Zhan g & Huil 120 05) showed 
how to derive the mass function corresponding to such barriers us- 
ing a simple numerical scheme. But for our purposes the Press- 
Schechter form suffice s. The most imp ortant prop erty is the density 
dependence. iBarkana & Loebl 120041) and iFurlanetto et alJ (2005) 
found that it is extremely similar to that contained in equation 
at least at high redshifts, and the former also showed that it fits 
simulations fairly well. 

Equations J4) and |5} give the dark matter halo abundance, 
which is not the same as the galaxy abundance. To connect the 
two we require the mean number of galaxies per halo, (N(mh)), 
the first moment of the halo occupation distribution. For simplicity, 
we will use the universal form of lKravtsov et all 12004). They split 
the mean occupation number into two parts. The first, (JV c (m/,)), 
represents the number of central galaxies in the halo and is unity 
if rrih > m m i n and zero otherwise. We will normally choose 
mmin by comparison to some observational detection threshold; 
note therefore that it need not reflect the actual minimum mass of a 
galaxy. The second part, {N s (mh)), describes the number of satel- 



lite galaxies above the same detection limit. iKravtsov et ail 120041) 
found that this has a simple form 



{N s (m h )) 



( m h V 

V C m m in / 



(6) 



From fits to the subhalo distribution in TV-body simulations, they 
estimate that C ~ 30 and (3 m 1 at z = over a broad range of 
halo masses. We will use these fiducial values in most of our calcu- 
lations. In order to estimate the importance of the halo occupation 
distribution, we will also compare to a model with {N(mh)} = 1. 

The total comoving number density of (observable) galaxies 
within a region is therefore 



n c g (m nlin \S v , M v ) 



dm h (N(m h )) n h (m h \S v , M v ), (7) 



where the superscript "c" denotes a cumulative quantity. We wish 
to compare to voids selected through galaxy surveys. We will there- 
fore set m m i n by comparison to the galaxy number densities in the 
Sloan Digital Sky Survey 1 and the 2dF Galaxy Redshift Survey. 2 
Specifically, we match to the number density of galaxies with r- 
band luminosity great er than a sp ecified value, using the luminosity 
function of Blan ton et all J2003). and to the c orresponding quanti- 
ties from the &j-band luminosity function of ICroton et all 120051) . 
The appropriate mass thresholds m m i n for the Press-Schechter and 
Sheth-Tormen mass functions are shown in Table Q We will also 
show some results for m m i n = 10 10 Mq in order to mimic an 
exceptionally deep survey. Note that normalizing in this way - to 
the cumulative number density of galaxies - tends to wash out dif- 
ferences between mass functions, insulating us from much of the 
uncertainty in using the Press-Schechter mass function. 



3.2 The linearized underdensity of voids 

The total (observed) galaxy underdensity in a void with physical 



size R v is 



1 + 5gal("lmin, 8 V , R v ) = 



n c a (m min \8^ , M v ) 
77 s ng(m min ) 



(8) 



Given m m i n , we solve this equation to find the linearized under- 
density required to produce a void of size R v and mean observed 
galaxy underdensity 5 ga i. 

Figure|2|shows the results of this inversion. The solid curves in 
panel (a) give the required underdensity for several different survey 
depths: M r < -20, -18, -16, and m h > 10 10 M , from top to 



1 See http://www.sdss.org/ 

See http://www.mso.anu. edu.au/2dFGRS/ 
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Figure 2. Linearized underdensity required to produce a mean (phys- 
ical) galaxy underdensity 5 ga j. (a): All curves assume <5 ga i = —0.8. The 
solid curves use our fiducial model, with M r < —20, —18, —16, and 
mj, > 10 10 Mq, from top to bottom. The dashed curves assume one 
galaxy per dark matter halo; we show M r < — 18 and mj, > 10 10 Mq. 
(b): The dot-dashed, solid, and dashed curves assume <5 ga i = —0.6, —0.8, 
and —0.9, respectively. We show results for M r < —20, —16 (top and 
bottom sets of curves). 



bottom. (Here and throughout, we will suppress the "—5 log h" in 
measured magnitudes for the sake of brevity.) We take 5 ga j = —0.8 
as a fiducial value. We plot 8^ as a function of mass variance a 2 . 
For reference, o 2 = 0.4, 1.9, and 4 correspond to comoving sizes 
R w 15, 5, and 2.5ft -1 Mpc, respectively. Note that this is not 
the physical size of the void (which is larger by a factor rj and 
varies with 5%), We halt each curve at M v = 2m m i n . The re- 
quired underdensity becomes more extreme in deeper surveys; this 
is because underdense regions are much more deficient in massive 
galaxies than small haloes in the extended Press-Schechter formal- 
ism. Haloes with masses mi < M v cannot collapse because too 
few density modes are available between (M v , mh) to pass the 
collapse threshold 5 C , but if a 2 3> o 2 (for small haloes) the loss of 
large-scale power becomes unimportant. In other words, the bias in- 
creases with galaxy mass, so the large galaxie s are more sensitive to 
the underlying density field lEfstathiou et al. 1988l lCole & KaiseJ 
ll989tlMo & WhitJl99d) . This is why 8% curves upward as a 2 in- 
creases, especially in the M r < — 20 case. There is also a slight 
trend for 8 V to increase as a — > 0; this is because the shape of the 
mass function changes in underdense regions, with more mass in 
small haloes. As the void size approaches infinity, slightly more of 
the mass shifts into galaxy-mass haloes. 

The two dashed curves in Figure |2j; show 8^ if we assume 
(N(m h )) = 1 for M r < -18 and m h > 10 10 Mq. This de- 
creases the threshold, because it increases the weight of the small- 
est haloes (which are least sensitive to the density). The value of 8^ 
becomes more sensitive to the halo occupation number for deeper 
surveys; a larger minimum mass pushes the survey to the steep part 
of the mass function, where the extra satellites in extremely large 
(but rare) haloes make only a small difference. Note as well that 



in this case the mh > 10 Mq threshold lies significantly below 
shell-crossing, where our model breaks down. 

Figure^ varies <5 ga i for M r < —20 and M r < —16 (top 
and bottom sets of curves). We let 5 ga i = —0.6, —0.8, and —0.9 
for the dot-dashed, solid, and dashed curves. This obviously has a 
strong effect on 8 V , because the massive galaxies are so sensitive 
to the underlying dark matter density. Note that the differences are 
largest at small a 2 (or large physical scales), because at smaller 
scales the finite size of the void limits the maximum halo mass. 
This is also why the differences are larger for the deeper survey. 

3.3 Haloes in voids 

Now that we have defined voids, we can look more closely at their 
resident halo populations. Figure [2] clearly shows that the intrinsic 
properties of voids depend sensitively on the mass limit of the sur- 
vey used to detect them: deeper surveys require larger dark matter 
underdensities and hence will presumably find smaller voids. But, 
in a given void, will smaller haloes be more abundant relative to 
the (brighter) galaxy population used to define the void? To answer 
this question, we let 



l + 8 h (m) = 



n h {m\8^,M v ) 
?7 3 n h (m) 



(9) 



We emphasize that Sh is the halo underdensity at mass m, not the 
galaxy underdensity. The analogous galaxy underdensity at a par- 
ticular mass is more difficult to calculate, because it requires the 
distribution of galaxy masses within each halo instead of simply the 
number of galaxies (i.e., how many satellite galaxies of mass m are 
contained in any galaxy group). Such models require considerably 
more machinery than is appropriate for our simple analytic model. 
We will therefore content ourselves with computing the halo dis- 
tribution, which will illustrate the most important features anyway. 
We will comment on its relation to the galaxy population below. 

Figure|3|shows predictions for the halo underdensity in voids 
with R v — 21 h^ 1 Mpc. The parameters in the two panels are iden- 
tical to those in Figure [2] In panel (a), we see that small haloes 
become less underdense if void selection is performed with larger 
galaxies. Within any given void, we also see that small haloes are 
less underdense than large haloes. This is, of course, because larger 
haloes are more biased than small ones. The disparity can be quite 
large for the M r < — 20 sample, with small haloes only half as un- 
derdense as those with rrih ~ 10 13 M© . The dashed curves show 
that these conclusions are qualitatively independent of (JV(mj)). 
In panel (b), we see that varying 5 ga i also modifies the fractional 
overdensity. However, its effects on the shape of the curve are rela- 
tively modest except at the most massive end. 

The horizontal bars on the right axis of Figure [5] show the 
physical dark matter underdensity within each void. Interestingly, 
it takes a relatively narrow range of values in panel (a), clustered 
around 8 ~ —(0.7-0.8), with \8\ decreasing as m m i n increases. 
The halo populations (especially at large masses) spread over a 
much broader range in Sh- This is essentially a consequence of 
the shape of 8(8 L ) shown in Figure Q a narrow range in physi- 
cal underdensities spans a large range of 8 L , which is the quantity 
relevant to our model for halo abundances. In general, 8 ~ 8h for 
haloes selected near m m i n . The density of haloes with m <C m m i n 
is closer to the mean because their bias is less than unity. 

Our model also predicts that Sh at a fixed mass will vary with 
the size of the void, because 5^ is also a function of scale. We show 
the implications in Figure|4j in this case we consider surveys with 
M r < —19 and M r < —16 (top and bottom sets, respectively). 
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Figure 3. Underdensity of haloes within voids of physical size R v = 
21 h^ 1 Mpc. (a): Assumes 5 ga i = —0.8. The solid curves use our fidu- 
cial model, with M r < -20, -18, -16, and m h > 10 10 M , from top 
to bottom. The horizontal bars on the right show 8 in each of these voids 
(same ordering from top to bottom). The dashed curves assume one galaxy 
per dark matter halo; we show Mr < —20 and mi > 10 10 Mq. (b): 
The dot-dashed, solid, and dashed curves assume 5 ga i = —0.6, —0.8, and 
—0.9, respectively. We show results for M r < —20, —16 (top and bot- 
tom sets of curves). The horizontal bars show the corresponding S for the 
M r < -20 cases. 
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Figure 4. Underdensity of haloes within voids of <5 ga i = —0.8 and several 
different radii. The upper thick and lower thin sets of curves assume M r < 
— 19 and —16, respectively. The horizontal bars at left show S for R v = 
7 h~ 1 Mpc voids. 



For large voids, R v > 7 h^ 1 Mpc, we see only a slight steepening 
of the curves, which occurs because the finite mass of the void lim- 
its the maximum halo mass. Thus, galaxy populations in voids with 
R v > 7 h" 1 Mpc should be nearly independent of the void radius. 
This is a consequence of the flatness of 5 V at small a 2 in Figure[2] 
which occurs because so little power exists on such large physi- 
cal scales. However, when R v ~ 2.1 h" 1 Mpc, the shapes change 
dramatically because the finite size of the void strongly limits the 
haloes inside. A comoving volume with radius R = 2.1 h^ 1 Mpc 
contains ~ 10 13 Mq, which is only ~ 20-100 times larger than 
the mass limits of these surveys. Thus the region requires only a 
modest underdensity to decrease the abundance of massive galax- 
ies. Such "voids" therefore sit near the mean density and smaller 
haloes (for which finite size effects can be neglected) will have 
nearly their mean density. This simplified model suggests that many 
properties of voids will depend on their spatial extents, at least if 
R v <7h- X Mpc. 

One consistent theme of this section is that large haloes are rel- 
atively less abundant than small haloes within voids. We can com- 
pare this expectation to observations, which suggest that the mean 
number density n ga i and the characteristic luminosity L* decrease 
significantly in voids but that the faint-end slope of the luminosity 
function remains nearly constant iCroton et al]|2005t iHovle et alj 
2005). Unfortunately, the comparison is not trivial. These studies 
measured the galaxy environment through the mean density within 
spheres of size ~ 8 h^ 1 Mpc (generally by comparing to the den- 
sity field of a volume-limited sample of galaxies). This is just large 
enough that finite size effects should not be dramatic, but it is still 
in a regime in which Poisson fluctuations in the galaxy distribution 



cannot be neglected. As a result, some of the "low-density" regions 
may in fact be at the mean density, which would wash out the void 
effect. An ideal test would focus on those galaxies identified to be 
within large voids (for example using the Voidfinder algorithm). 
Moreover, we must also bear in mind that Sh ^ <5 ga i. In particular, 
many small galaxies will be satellites inside large haloes; such a 
population will be efficiently suppressed in just the same way as 
bright central galaxies. Thus we expect the galaxy luminosity func- 
tion in voids to be flatter than the halo mass function. 

With these caveats, our model naturally explains the decrease 
in n ga i and L*, especially because these regions are relatively small 
and finite size effects help suppress galaxies with m > 10 12 Mq. 
Note that this probably explains the discrepancy between the ob- 
served Li, and that predicted by Coorav 1 2005), who let R v — ► oo 
rather than including a cutoff at R v ~ 8/i _1 Mpc which would 
have quenched the formation of massive galaxies. On the other 
hand, our models robustly predict an excess (by about a factor of 
two) between the abundance of small and large haloes in voids. 
If there were a one-to-one correspondence between halo mass and 
galaxy luminosity, and if the luminosity function could be approxi- 
mated by a power law n oc L~ a over this interval, we would expect 
a to be ~ 0.15 larger in voids than in mean density environments. 
The suppression of satellite galaxies would reduce the di sparity. 
Even so, our result is within the error bars of lCroton et alJ f2005), 
though their best-fit values indicate no significant change. Thus, 
there is currently only weak evidence, at best, for non-standard 
galaxy formation within voids. CDM models robustly predict an 
excess of small haloes inside of voids, but it is not particularly 
strong and requires care to interpret. 
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4 VOID ABUNDANCES 



4.1 The excursion set abundance 



Sheth & van de Weveaertl <2004 used the excursion set formalism 
to estimate the number density of voids. They defined a void as a 
region with dark matter underdensity 5% ■ At first sight, the void 
abundance seems to follow from equation 0, but with 8^ — > 
8^: this is the same problem as the halo abundance, except that 
the absorbing barrier is negative rather than positive. However, 
ISheth & van de Weveaertl 12004) pointed out one crucial differ- 
ence, which they called the "void-in-cloud" problem. Consider a 
small region with 8 L < 8^ that is contained in a larger region with 
8 L > 8c . Equation would have assigned such a point to a small 
void; however, physically we know that this "void" lies inside of 
a collapsed object and hence has been crushed out of existence. 
(This problem does not occur for nh because haloes are allowed to 
collapse inside voids.) 

Thus, the appropriate diffusion problem has two barriers: we 
wish to compute the first-crossing distribution for the void barrier 
8^ while excluding trajectories that have already crossed a barrier 
with 8p > that describes void-crushing. The simplest assumption 
is to take both these barriers to b e independent of a 2 ; in that case 
ISheth & van de Weveaertl j2004) showed that the mass function is 
(see also the alternate derivation in the Appendix) 



i v (m v ) = — j 



d In a 



dlnm„ 



E 



x exp 



nW 

~i{8\iioy 



u 2 tt 2 D 2 sm(mvD) 
(8*; /a) 2 nn 



where 



D = 



(10) 



(11) 



describes the relative importance of void-crushing; a fraction (1 — 
D) of all matter lies inside voids in this model. Note that 8^ and 8p 
implicitly depend on redshift, because we evaluate a at the present 
day. 

The lshefh & van de Weveaertl h004 model contains one cru- 
cial assumption - that voids can be defined at a constant 8^, in- 
dependent of scale - and two free parameters. In what follows 
we will usually set the void-crushing parameter 8^ — 1.06, 
which is the linearized overdensity at turnaround according to the 
spherical collapse model. This is reasonable because it marks the 
point at which larger scale overdensities will begin to collapse 
around their resident voids, but we will also let 8^ take larger 
values for illustrative purposes. The second free parameter is 8„; 
ISheth & van de WeveaerT<2004 argued that 5% = -2.8 (shell- 
crossing) is an appropriate choice, because that marks the maxi- 
mum point of efficient expansion. 

However, we have already seen that the dark matter underden- 
sity is generally significantly deeper than the galaxy underdensity. 
Instead, if one wishes to compare to voids observed in redshift sur- 
veys, 8^ should be chosen from the 5 ga i and m m i n appropriate to 
a given void search. The relevant absorbing barriers are precisely 
the curves shown in Figure [2] only in voids selected from excep- 
tionally deep galaxy surveys do we expect shell-crossing to have 
occurred. Interestingly, the barriers are also relatively flat, at least 
for reasonably deep surveys and on large scales. We will therefore 
take 8v = constant and use equation I lOi to compute the number 
density of voids. We evaluate 8^ at R v = 21 h" 1 Mpc in the fol- 




Figure 5. Abundance of voids in the excursion set formalism. All curves 
follow our fiducial model; they vary the minimum detectable galaxy mass 
as shown in the legend. Panels a and b show differential and cumulative 
filling factors of voids, respectively. 



lowing; we examine how n v (m v ) varies with the barrier shape in 
j44lbelow. 

Figure [5] shows the void abundance for several different sur- 
vey depths. Panel a shows the fraction of volume filled by voids 
of a given radius; note that we use the physical void radius here, 
including the extra gravitational expansion. Panel b shows the vol- 
ume filling fraction of voids, 



F vo i(> R) = rf 



am v — n v 



(12) 



where the lower integration limit is the mass corresponding to a 
void of physical radius R; the prefactor rf accounts for the extra 
gravitational expansion. We halt each curve when the void mass is 
less than twice m m i n . 

Figure [5] makes clear several key characteristics of the void 
size distribution. First, n v has both large and small mass cutoffs. 
The former occurs where a ~ 8^ , above which voids are ex- 
ponentially suppressed because of the lack of large-scale power 
(the same mechanism that provides the large-scale cutoff in the 
halo mass function). The small mass cutoff occurs because of the 
void-crushing threshold 8^ . The resulting characteristic scale R c 
is quite sensitive to the galaxies used to select voids: it ranges 
from R c w bh' 1 Mpc for M r < -16 to R c » Wh' 1 Mpc for 
M r < — 21, nearly an order of magnitude in volume. The shape, on 
the other hand, does not change much across the different selection 
thresholds. Finally, we also see that F vo i > 1: the voids appar- 
ently fill a volume larger than the universe. This occurs because we 
have let each and every void expand to its full size, even those sur- 
rounded by relatively weak overdensities just short of turnaround. 
These filling fractions should thus not be taken too seriously; they 
suggest only that, when voids are selected from galaxy surveys, we 
expect empty regions with R^Sh^ 1 Mpc to fill a large fraction of 
the universe. 

The largest possible voids also provide an interesting bench- 
mark for comparison to surveys l Blumenthal e t alJll992t) . To es- 
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timate J? max , we find where n v = (c/Ho)~ 3 ; i.e., the size for 
which we expect one void per Hubble volume. This will likely over- 
estimate the maximum size in a real survey, because voids were 
smaller in the past (see i|4,3l below). Of course the maximum size 
will also vary with the magnitude limit of the survey; for M r < 
(-16, -19, -21) we find i? max « (28, 32, 45) ft -1 Mpc. These 
are compatible with existing redshift sur veys; the largest voids i n 
the 2dF survey have R ~ 25/i _1 Mpc faovle & Vogelevll20o4) . 
We will make a more detailed comparison to observations in i |4.5l 
below. 

Finally, we note some simple properties of F vo \. Ignoring 
void-crushing, the void filling fraction becomes 



CO R) 



rf erfc 



V2a{R) 



(13) 



This is the same expression found by Dubi nski et al] fl993). except 
that they assumed Sy corresponded to shell crossing and included 
an extra factor of one-half (which occurred because they ignored 
those trajectories that joined larger voids). Void-crushing is unim- 
portant for voids larger than the characteristic size (see 34.4l below). 
so equation 1 1 3i provides a reasonable estimate for F vo i and shows 
explicitly how it depends on 5„ (and hence on the survey charac- 
teristics). Dubin ski et alJ fl9 93) pointed out the interesting coinci- 
dence that _F vo i w 1 at the present day; this remains true in our 
model. Only recently have large voids in the faint galaxy popula- 
tion come to dominate the geography of the universe. 



(a) Thick: M r <-16 
Thin: M <-20 




H-+ 



10= 



Figure 6. As Fig. [5] except we vary the mean galaxy underdensity within 
the voids: 5 ga i = —0.8, —0.9, and —0.6 (solid, dashed, and dot-dashed 
curves, respectively). The thick and thin sets of curves take M r < —16 
and M r < -20. 



4.2 Void depth and abundance 

We have seen that one key quantity in defining voids is the mean en- 
closed galaxy underdensity <5 ga i . Figure[6]shows how the void sizes 
vary with this criterion. The curves correspond to the same cases 
shown in Figure |2l?. We see that <5 ga i has an enormous effect on the 
characteristic void size R c : it changes by a factor of three over this 
range of galaxy density (or a factor of nearly thirty in volume!). 
Moreover, in panel (b), we see that large, deep voids (<5 ga i = —0.9) 
only fill a small fraction of space. 



4.3 Voids at higher redshifts 

Just as with the halo mass function, the excursion set formalism 
also allows us to study how n v (to) evolves with redshift. We show 
the predicted void distribution at z = 1 in Figure^ In order to fa- 
cilitate comparison with the z — results, we hold m m in constant 
between z = 1 and the present day (as in Table 1). Of course, be- 
cause the galaxy luminosity function evolves over this interval, the 
number density of these haloes evolves over this redshift interval. 

The voids also evolve, with the characteristic radius increas- 
ing significantly (by about a factor of two in each case) and the vol- 
ume filling fraction also increasing (again by a factor of ~ 2-3). 
For example, only ~ 10% of the universe lies inside of voids with 
R > 10 h" 1 Mpc selected from m h > 1.5 x 10 12 M Q galaxies 
(as opposed to ~ 33% at z = 0). This implies that voids con- 
tinue to grow, albeit relatively slowly, even after the cosmological 
constant dominates the energy density. O ur results are in qualita- 
tive agreement with lConrov et alJ j2005fl . who find that voids are 
both smaller and rarer at z ~ 1. A precise comparison is difficult 
because they use statistical techniques (specifically the void proba- 
bility function) rather than identifying individual voids. 




R (h" 1 Mpc) 

Figure 7. As Fig.|5] except for 2 = 1. The galaxies are selected with the 
same mass thresholds as in Fig.|5]or Table Q] (with the solid curve having 
*rtmin > 10 10 Mq and mass increasing to the right). Note that the z = 1 
number densities differ from the z = values. 

4.4 Model assumptions 

Unfortunately, this model for n v (m) contains a number of free pa- 
rameters and simplifications, so it is not nearly as well-specified as 
the halo mass function nh(jn). Here we will examine how sensitive 
our results are to these assumptions. We begin with the halo occu- 
pation num ber (N(rrih))- The sol id and dashed curves in Figure|8}i 
contrast the lKravtsov et alj (;2004) value (our fiducial model) and a 
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Figure 8. Void abundance for several parameter choices. All curves take 
our fiducial model except as noted, (a): Dashed curves assume one galaxy 
per dark matter halo, while the dotted curves use the Sheth-Tormen mass 
function (b): The solid, dashed, and dotted curves set the "void-crushing" 
parameter Sp = 1.06, 1.69, and 10, respectively. In both panels, the three 
sets of curves assume surveys with m m i n = 10 10 Mq, M r < —18, and 
M r < -20, from left to right. 



Figure 9. Void abundance for different barrier prescriptions. All curves take 
our fiducial model except as noted, (a): Solid, dashed, and dotted curves set 
the (constant) excursion set barrier by evaluating S v at R v = 70, 7, and 
2.1h — 1 Mpc, respectively. (The fiducial model uses R v = 21h~ 1 Mpc.) 
(b): The dashed curves use a linear fit to the excursion set banier to compute 
the void abundance. In both panels, the three sets of curves assume surveys 
with M r < -16, -18, and -20, from left to right. 



model with one galaxy per halo. We have previously seen that (TV) 
only significantly affects 8y if m m i n is relatively small; as a result, 
n v (m) is only substantially affected for the m m i n > 10 10 Mq 
survey. In this case the voids shrink by about 40% in radius, with 
R c ~ 3ft" 1 Mpc. Thus, if we selected voids from the halo un- 
derdensity (as, for example, may be possible by associating groups 
of galaxies with single haloes), we would expect slightly smaller 
voids. The difference is minimal for bright galaxies because the 
mass function is steep in that regime. 

Figure [SJ; also compares the void sizes computed from the 
Press-Schechter mass function with those from the Sheth-Tormen 
mass function. To perform this comparison, we have scaled the 
Sheth-Tormen mass function with de nsity in the same way as the 
Press-Schechter mass function (Barkana & Loeb 2004). Thus the 
only real difference is the mass threshold m m i n , although that dif- 
ference can be substantial. We see that n v (m) remains nearly un- 
changed. Thus while the void distribution may in detail depend on 
the exact form of (m), it is not likely to be substantially affected. 

Another unknown parameter in the model is 5~, 
whic h describes void-crushing . Figure [8^ shows its effects. 
ISheth & van de Weveaertl j2004) argue that Sp is unlikely to be 
smaller than 1.06, because voids should continue to increase in 
size until turnaround. We therefore consider Sp = 1.69 (where 
voids are only destroyed when they lie inside of collapsed haloes) 
and Sp — 10 (where void crushing is negligible). This parameter 
clearly does have an effect on n„(m): ignoring void crushing 
dramatically increases the number of small voids. However, it has 
no effect on large scales, because these voids cannot be inside of 
extremely dense regions anyway. It also has little effect on the 
characteristic size R c . Thus, while 5 p may affect the distribution 
of small voids, it is unlikely to affect the much more striking 
population of large voids. 



Perhaps the most important simplification we have made is 
to ignore the scale dependence of <5.„ so that we could use equation 
J 101 . Figure[2|shows that the threshold actually varies with the mass 
of the void. In principle, we should find the exact size distribution 
by following a procedure similar to lZhang & Huil 120051) . but with 
two absorbing barriers. However, it is relatively easy to see that in 
most cases the constant barrier approximation suffices for our pur- 
poses, though more sophisticated comparisons in the future may 
require a more exact solution. One straightforward test is to take 
5v to be constant but to evaluate its amplitude at several different 
scales. In Figures[3]|8] we have used S^(R = 2lh~ 1 Mpc). In Fig- 
ure^, we show the corresponding distributions for the threshold 
evaluated at R = 70, 7, and 2.1h~ 1 Mpc (solid, dashed, and dot- 
ted curves). We show results for M r < —16, —18, and —20 (left, 
middle, and right sets). Clearly, for deep surveys n v (m) is nearly 
independent of the scale at which we evaluate the barrier. This is, 
of course, because the corresponding thresholds in Figure [2jz are 
nearly independent of a 2 . 

However, if M r < —20, the barrier rises rapidly at small 
scales. From Figure the distribution is still nearly indepen- 
dent of scale so long as R v >7/i _1 Mpc, but using Sy(R — 
2.1h~ l Mpc) gives entirely different results. Of course, this choice 
is not self-consistent, because the distribution is dominated by 
much larger voids. To better estimate the effects of this "moving" 
barrier, we note that it is approximately a linear function in a 2 . We 
would therefore ideally like to solve for the mass function in the 
presence of a constant barrier Sp and a linearly increasing barrier 
5 V + /3a 2 . Unfortunately, we were unable to find an analytic solu- 
tion for that distribution (though a straightforward numerical solu- 
tion may be possible by extending Zhana & Hui 2005). We there- 
fore instead turn to a model in which the void-crushing barrier is 
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also linear with the same slope as the void barrier. In the Appendix, 
we show that the appropriate mass function is 

n„(m„) = n v (m v ) exp(—(3S^ — /3 2 a 2 /2). (14) 

(Note that this differ s from the solution presented by 
Sheth & van de Wevsaert 2004, which corresponds to a case 
with two barriers of opposite slopes; see the Appendix for a more 
detailed discussion.) Treating 5p as a linear barrier is probably 
not a bad approximation. For one, the amplitude of 5p does not 
strongly affect n v (m), as shown by Figure [8j>. Moreover, most 
trajectories that cross 8y would have crossed the void-crushing 
barrier at small o 2 , when /3a 2 <C 1. Finally, the appropriate 
choice of 5p is somewhat arbitrary anyway, so a linear model may 
actually be more accurate. 

The similarity of n v (m v ) and n v {m) suggests that the linear 
barrier will have only a small effect on the distribution. In the cases 
of interest, (3 > 0, so the first exponential factor will be larger than 
unity and more mass will lie inside of voids. This is intuitively ob- 
vious, because an increasing negative barrier is easier to cross. The 
second term affects the shape of the mass function, but it is only 
important for large f3 2 o 2 ; in the regime we study, this factor is al- 
ways smaller than unity. Figure[9j> confirms these expectations. For 
M r < —20, the shape remains nearly invariant (except at the small- 
est scales) but more mass lies inside of voids. Most importantly, the 
characteristic scale R c is nearly identical in the two cases. The lin- 
ear barrier makes no difference to fainter galaxy samples, because 
the corresponding barrier is so nearly flat. 

4.5 Comparison to observations 

We are now in a position to com pare our model to the ob- 
served distribution of voids. Hovle & Vogelev 1 2004) compiled 
the most complete sample to date from the 2dF Galaxy Redshift 
Surv ey. They use d the Voidfinder algorithm <El-Ad&Pirarll997t 
IHovle & Vogelevll2002l) to identify voids with R > 10 h' 1 Mpc 
in a volume-limited sample of galaxies with Mhj < — 19. s They 
found that ~ 35% of the universe is contained inside voids with 
R > Wh,- 1 Mpc; the mean effective radius of their sample was 
R w 15 h' 1 Mpc, with a maximum of R « 25ft -1 Mpc. Their 
algorithm separates galaxies into "wall" and "void" populations; 
the latter can sit inside of voids but the former are excluded. They 
measured a mean galaxy underdens ity 5 ea i = —0.93 w ithin voids. 

With this magnitude limit, the H ovle & Voselevl l2004h sam- 
ple should find voids somewhere between the two dot-dashed 
curves in Figure [5] or near the thin curves in Figure |6| Given the 
measured <5 ga j = —0.93, the most naive comparison is to the 5 ga i = 
—0.9 curve in the latter, which has F vo i (R > Wh' 1 Mpc) « 16% 
and R c ~ 5 hT 1 Mpc. Both of these are considerably smaller than 
the observed value, but it is unclear how significant the discrep- 
ancy is. Unfortunately, identifying voids at the predicted character- 
istic scale is difficult, because Poisson fluctuations in the galaxy 
distribution are substantial on Mpc scales. As a result, it is im- 
possible to know whether an observed galaxy underdensity with 
R ~ 5/i _1 Mpc corresponds to a true dark matter underdensity or 
a stati stical fluctuation. Largely for this reason, Hovle & Voeelev 
(2004) restrict their search to voids with R > 10 /i -1 Mpc. We 
therefore do not necessarily consider the small R c predicted by our 

3 IHovle & Vogeiel (2002) do not explicitly state the absolute magnitude 
limit for their main sample, but comparison of its size to their Table 2 yields 
this value. 



model to conflict with observations; however, the discrepancy in 
Fvoi f° r large voids is more worrying. We do remark that a large 
population of small voids have clear observational signatures, in- 
creasing the number of ~ 5k' 1 Mpc empty or nearly empty re- 
gions well above those expected in a random distribution. This 
would manifest itself in such statistical measures as the void prob- 
ability function (the probability that a region has zero galaxies; 
IWhitd l 1979) or the underdensity probability function (the proba- 
bility that a region has a given underdensity; IVogelev et all 1989), 
but such statistics are (at least for the moment) beyond the capabil- 
ities of our model. 

There are several more subtleties in a comparison to obser- 
vations. Most importantly, the Voidfinder algorithm does not use 
a density threshold, so the appropriate choice of <5 ga i is not obvi- 
ous. The only "void" galaxies allowed in this approach are those 
that are well-isolated. But not all galaxies inside real voids need 
be isolated: if some fraction of the wall galaxies are actually con- 
tained inside regions with 5 L < S„, the true voids would be less 
underdense than allowed by the algorithm. Indeed, the <5 ga j = — 0.8 
curves in Figure HJpredict F mi (R > 10 h' 1 Mpc) « 25-33% and 
R c ~ 7-8 hT 1 Mpc, much closer to the observed values. There are 
solid physical reasons to expect confusion between genuine void 
and wall galaxies. Spherical underdensities tend to evolve by evacu- 
ating their in teriors, wit h mass accumulating in a shell near the edge 
(see, e.g- lDubinski et alJl993h . lGotflober et alJl2003l) showed that 
massive haloes pile up at the edges even more strongly. Thus it is 
reasonable to expect that many of the bright galaxies belonging to 
a large-scale underdensity lie near its edges, where they are diffi- 
cult to identify unambiguously. Moreover, real surveys operate in 
redshift space. Because voids are expanding more rapidly than the 
Hubble flow, their observed volumes will be larger than the physi- 
cal volumes used by our model. Gol dberg & Vogelevl 120041) show 
that the enhancement can be > 20% near shell-crossing. 

An equally important consideration is the geometrical method 
used by Voidfinder to identify voids. It fills the gaps between wall 
galaxies with spheres and then merges overlapping spheres into 
discrete entities. This procedure may join two neighboring voids 
that our approach would consider distinct, especially if the wall 
between them is weak or if we include Poisson fluctuations. Galax- 
ies on void walls tend to es cape along the walls when neighbor- 
ing voids merge toubinskietafl ll993). making this picture even 
more likely. Finally, voids are not uniform; rather, they are com- 
posed of a patchwork of (deeper) voids separated by weak filaments 
(Gottlober et al. 2003). In the excursion set picture, this is the ana- 
log of the progenitor distribution of collapsed objec ts and could in 
princ iple be computed through similar techniques lLacev & Cold 
1993), with the important difference that "sub- voids" can more 
easily retain their distinct identity as deeper underdensities in the 
galaxy distribution. The interaction of this structure with void se- 
lection algorithms is not trivial. 

Another clear predic tion of our model is tha t voids become 
larger as m m i n increases. Hovle & Vogelev (2004) addressed this 
question with statistical techniques. Specifically, they computed 
the void probability function and underdensity probability function 
for volume-limited samples with a range of magnitude thresholds 
(fromM 6 / < -16toM{,j < -21). They found that both the frac- 
tion of the universe filled by voids and their characteristic size in- 
crease with the luminosity threshold, in qualitative agreement with 
our results. A more precise comparison requires an extension of our 
model to these statistics. 

In summary, although our model predicts a significant popu- 
lation of large voids, they are still somewhat smaller than the ob- 
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served structures. This could be a real discrepancy - in which case 
voids clearly require more sophisticated modeling - or it could be a 
combination of the contrasting selection techniqes, the internal void 
structure, and redshift-space effects. However, our qualitative con- 
clusions (large voids, increasing in size as the underlying galaxy 
population brightens) do match the observed trends, a significant 
step forward for analytic models. Unfortunately, precise compar- 
isons to observations will still require numerical simulations. 

4.6 Comparison to simulations 

Another useful co mparison is to voids i nside cosmological simula- 
tions. In particular, B enson et alj 1200 3l) used the Voidfinder algo- 
rithm to identify voids inside an iV-body simulation. They used a 
semi-analytic model to place galaxies within dark matter haloes and 
selected the voids based on the predicted properties of the galaxies. 
They presented the void size distributions for several different sur- 
vey depths. Our results compare favorably with theirs. Like us, they 
found larger voids for brighter galaxies as well as larger voids for 
galaxies compared to the dark matter (though below ~ 10 h^ 1 Mpc 
the two appear to converge; this may be related to the increasing 
importance of Poisson fluctuations). They also found that the size 
distribution only peaks at R > 10 hT 1 Mpc for the brightest galax- 
ies {M r < —21.5 and Mbj < —20.5 in their models); for fainter 
galaxies their distributions continue rising to R ~ 7 Mpc (be- 
low which their selection technique becomes incomplete). Finally, 
they find that voids typically have 8 w —0.8 and <5 ga i ~ —0.9, 
with both quantities decreasing rapidly toward the center of the 
void. The density increases sharply near, but slightly beyond, the 
nominal void radius. This provides a hint that the density structure 
may play a role in the Voidfinder selection, making direct compar- 
isons with our predictions somewhat more difficult. Nevertheless, 
the good qualitative, and even reasonable quantitative, agreement 
of our model with these simulations provides strong support for the 
maj or features of our ap proach . 

IColberg et alj 120051) recently examined voids in a series of 
large-volume iV-body simulations. They found voids to be much 
smaller than our predictions, with 90% of the volume contained 
in voids with R < 2.5 h^ 1 Mpc. However, they selected voids 
based on a dark- matter density criterion 8 V < —0.8 (corresponding 
to shell-crossing). This is much more restrictive than our galaxy- 
density criterion. The best comparison is to the mi > 10 10 Mq 
survey in Figure|3| whose 8% is close to shell-crossing (see Fig.|5J. 
For this sample, we find R c w 4/i _1 Mpc, with much of the uni- 
verse in smaller voids. The remaining discrepancy may result from 
the search algorithm. They began their void searches around local 
minima in the density field, extending them until the enclosed den- 
sity exceeded 8 V . Such a prescription may tend to find "sub-void" 
progenitors rather than the larger structures that correspond to the 
observed voids in the galaxy distribution. 



5 DISCUSSION 

We have described a simple analytic mo del for voids in the galaxy 
distribution. Our model is based on ISheth & van de WeveaerJ 
(2004), who showed how to apply the excursion set formalism to 
underdensities. Its most important parameter is 8 V , the linearized 
underdensity of a void. Those authors originally set 8^ to be the 
density corresponding to shell-crossing. However, this condition 
produces voids much smaller than the observed structures, whose 
masses can exceed those of galaxy clusters lPiran etaJl99l . Our 



major contribution has been to show how to define 8 V through the 
galaxy underdensity. Because galaxies are biased relative to the 
dark matter, voids can be nearer to the mean density than the ob- 
servations naively indicate. This significantly increases the charac- 
teristic void size and the volume filling fraction. 

Our model predicts voids with characteristic radii R c ~ 5- 
10/i -1 Mpc. This is similar to, though somewhat smaller than, the 
voids found in galaxy redshift surveys iHovle & Vogele^l2002l) . It 
is a closer match to the void population in semi-analytic galaxy for- 
mation models iBenson et alj|2003h . Because bright galaxies are 
more highly biased, our model also predicts that voids selected 
from shallow surveys s hould be c haracteristically larger, in agree- 
ment with observations l Hovle & Vogelev 2004). However, we also 
predict larger galaxy densities inside voids than observed. The sig- 
nificance of these discrepancies is not clear: we argued in i|4.5l that 
the algorithms used to identify voids in surveys, redshift-space dis- 
tortions, and the internal structure of voids are all important in de- 
tailed comparisons to the observations. The last of these is espe- 
cially crucial, because simulations sh ow that galaxies insid e voids 
tend to congregate near their edges iGottlober et aljEo Q31. where 
they are difficult to separate from genuine "wall" galaxies. Two 
other subtleties may also affect the comparison. The first is the non- 
linear evolution of substructure within the void; we have used lin- 
ear theory to describe the halo population inside voids, which likely 
breaks down in some regimes because power can be transferred be- 
tween s cales. The second is Poisson var iation in t he galaxy n umber 
counts ISheth & LemsorJll99St ICasas-Miranda et alj [2002). This 
effect makes direct detection of small voids impossible, forcing ex- 
isting surveys to search only for voids with R > 10ft _1 Mpc. 

Thus it is difficult to make a precise comparison between our 
simple model and observed voids. However, the qualitiative agree- 
ment is reasonable. Our analytic approach is the first to contain 
voids with characteristic radii ~ 10 h~ Mpc and hence to be in 
even qualitative agreement with the observations. Our model shows 
that the enormous extent of observed voids is not particularly sur- 
prising and should not be perceived as a "crisis" for the CDM 
model. 

By using the excursion set approach, we have also self- 
consistently predicted the halo population inside of voids. In agree- 
ment with naive expectations, small haloes should be less under- 
dense than massive haloes. But the discrepancy is by no means 
large - typically smaller than a factor of two over three decades 
in halo mass. Including satellite galaxies in this calculation (which 
we have not done) will decrease the difference. Thus, the predicted 
steepening of the mass function in low-density environments is 
only modest and is not rul ed out by existing observations of th e 
galaxy luminosity function ICroton et all200 j ; lHovle et all20o3) . 
We also emphasize that these observations determined the environ- 
mental density on relatively small scales (~ Sh^ 1 Mpc), where 
random fluctuations in the galaxy field may mimic true underdensi- 
ties. If so, differences between mean and low-density environments 
will be further washed out. Moreover, on such scales, finite-size ef- 
fects become important in setting the characteristic luminosity L* . 
The test can be sharpened by focusing on the luminosity function 
within large, easily identified voids - although, even there, sepa- 
rating void and wall galaxies may be difficult. Thus, we find no 
compelling reason to believe that galaxy formation differs inside 
and outside of voids - although we have no evidence against such 
a possibility, either. 

Statistical properties of the galaxy distribution can also be 
used to test our predictions, especially for small voids that cannot 
be identified unambiguously. The void probability function quanti- 
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fies the probability that a sphere of a given size contains no galaxies 
IWhitel 19791) . and the underdensity probability function quantifies 
the probabi lity that a sphere is more underdense than some speci- 
fied value jVogelev et aIT"9 89). These measures are free from bias 
in the void selection process, but they are harder to predict directly 
from the excursion set formalism and are less useful in identifying 
galaxies that reside in voids. Our model is in qualitative agreement 
with the observed trends in these statistics l Hovle & Voaelev 2004 ; 
IConrov et ai]l20 05). but more work is needed to connect our for- 
malism to them. 
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APPENDIX A: THE FIRST CROSSING DISTRIBUTION 
FOR TWO BARRIERS 



ISheth & van de Weygaert! 120041) showed how to derive equation 
\\{)\ using the Lapla ce transform. Here w e present an alternate 
derivation (similar to McQ uinn et all2 005) and generalize it to the 
case of two linear barriers with equal slope. We begin by consider- 
ing two constant barriers Sp — B p and 8y — B v . The distribution 
of trajectories Q(S L , S), where S = a 2 , obeys the diffusion equa- 
tion 



dQ 
dS 



1 d 2 Q 



2 d(S L 



(Al) 



The appropriate boundary conditions are Q(B P ,S) = 
Q{B V ,S) = and Q(S L ,0) = 5 D (S L ), where S D is the 
Dirac delta function. We first make the simple transformation 
y — S L — B v . We then assume Q(y,S) = f(S)g(y); equation 
JA1> becomes 



f dS 



2g dy 2 



(A2) 



where a is a constant. The general solutions are f(S) = 
exp(— a 2 S) and g(y) = exp(— \f2ay). The boundary conditions 
g(0) = g{B p — B v ) = select the sine function for the latter and 
force its argument to take discrete values; thus 



Q(y, S) = a„ sin 



niry 
B v — B v 



exp 



n 2 n 2 S 



2(B P — By) 



.(A3) 



We fix the constants a n by matching to the function Q(y,0) 
S D {y + B v ). Thus 



Q(y,s) = 



2 sin(n7rD) 



n=l 



B D 



x exp 



— B v 

nVff 
2(B P — By) 



niry 

Bp — By 



(A4) 



where D = —B V /(B P — B v ). The total rate at which trajectories 
disappear from the permitted region is 



d_ 

dS 



dyQ(y,S) 



(i 



dy 



(A5) 



where we have used equation <AH . We interpret the two terms on 
the right hand side as the rate at which trajectories flow across the 
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two barriers. Thus the first-crossing distribution for the void barrier 
is 



F V (S) = 



dQ 
dy 



E 



(mrD) 2 sin(mvD) 



exp 



„2 2 n 2 

n it D 



2BI/S 



(A6) 



which yields equation <10i after converting to mass units. 

We now consider the case of two absorbing barriers linear in 
S, Sp = B p + SiS and 6% = B v + BiS 1 . We make the coordinate 
transformation u — Bi(8 L — BiS) — B\B V ; then the boundary 
conditions on u become Q(0, S) = and Q(ui, S) = 0, where 
mi = Bi(B p — B v ). The boundary condition on S becomes 

Q(u,0) = 5 D {5 L ) = B 1 5n{u + BxB v ) exp(-u - BiB v ).(AT) 

The diffusion equation IA1I becomes 



ds 



Bf 
2 



+ 2^ 
du 2 du 



(A8) 



We again assume a separable solution Q(u, S) = f(S)g(u). The 
general solutions are exponentials, and the boundary conditions on 
g select the sine function with discrete arguments. Thus the solution 
is again a Fourier sine series in u. The Fourier coefficients a n can 
be fixed by matching the solution to equation <A7t . We then find 
that 



Q(u,S) 



E 



2 sin(n7rD) 
B„ — B v 



x exp 



—B±B V — 



BlS 



B\ [Bp — B V/ 

2 2 c 

nub 



2(B P - B v ) 2 



;(A9) 



note the similarity to equation <A4> . The rate at which trajectories 
cross the barriers is 



dQ 



du 



BiIQlo 1 • (A10) 



But Q(0, S) — Q(ui, S) = 0; thus we can again identify the first- 
crossing distribution of the linear void barrier F* with the u = 
component of the first part, so 

Bx dQ 

o 



fUs) 



2 du 



F V (S) exp(—BiB v — BfS/2), 



(All) 



from which the mass function follows directly. Note that this ex- 
pressio n differs from equation (CIO) of ISheth & van de Wevgaern 
(2004). Those authors tried to derive the same distribution through 
Laplace transforms. Their solution required the first-crossing dis- 
tribution of Sp and 5% to be identical (their equation C6); this is 
not the case if the barriers have equal slope B\ , because then they 
are not symmetric about the 5 L — axis. Instead their solution ap- 
plies to a case in which the slopes of the two barriers have opposite 
signs. 
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